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We show that combined electrostatic and radiative fields can greatly amplify the directional 
properties, such as axis orientation and alignment, of symmetric top molecules. In our computational 
study, we consider all four symmetry combinations of the prolate and oblate inertia and polarizability 
tensors, as well as the collinear and perpendicular (or tilted) geometries of the two fields. In, 
respectively, the collinear or perpendicular fields, the oblate or prolate polarizability interaction 
due to the radiative field forces the permanent dipole into alignment with the static field. Two 
mechanisms are found to be responsible for the amplification of the molecules' orientation, which 
ensues once the static field is turned on: (a) permanent-dipole coupling of the opposite-parity 
tunneling doublets created by the oblate polarizability interaction in collinear static and radiative 
fields; (b) hybridization of the opposite parity states via the polarizability interaction and their 
coupling by the permanent dipole interaction to the collinear or perpendicular static field. In 
perpendicular fields, the oblate polarizability interaction, along with the loss of cylindrical symmetry, 
is found to preclude the wrong-way orientation, causing all states to become high-field seeking with 
respect to the static field. The adiabatic labels of the states in the tilted fields depend on the 
adiabatic path taken through the parameter space comprised of the permanent and induced-dipole 
interaction parameters and the tilt angle between the two field vectors. 



I. INTRODUCTION 

Directional states of molecules are at the core of all 
methods to manipulate molecular trajectories. This is 
because only in directional states are the body-fixed mul- 
tipole moments "available" in the laboratory frame where 
they can be acted upon by space-fixed fields. 

In the case of polar molecules, the body-fixed perma- 
nent dipole moment is put to such a full use in the lab- 
oratory by creating oriented states characterized by as 
complete a projection of the body-fixed dipole moment 
on the space-fixed axis as the uncertainty principle al- 
lows. 

The classic technique of hexapole focusing selects pro- 
cessing states of symmetric top molecules (or equivalent) 
which arc inherently oriented by the first-order Stark ef- 
fect [TJ , independent of the strength of the electric field 
applied (at low fields). For the lowest precessing state, 
only a half of the body-fixed dipole moment can be pro- 
jected on the space fixed axis and utilized to manipu- 
late the molecule. Pendular orientation is more versa- 
tile [2J applicable to linear molecules as well as to 
symmetric and asymmetric tops, but works well only for 
molecules with a large value of the ratio of the body- 
fixed electric dipole moment to the rotational constant. 
However, for favorable values of this ratio, on the order 
of 100 Debye/cm -1 , the ground pendular state of a : S 
molecules utilizes up to 90% percent of the body-fixed 
dipole moment in the space-fixed frame at feasable elec- 
tric field strength, of the order of 100 kV/cm. 

Pendular states of a different kind can be produced by 
the induced-dipole interaction of a nonresonant laser field 
with the anisotropic molecular polarizability |4j [5] [6j [7] . 
The resulting directional states exhibit alignment rather 
than orientation, and can be used to extend molecular 



spectroscopy |5], suppress rotational tumbling [91 ITUl ITT] . 
focus molecules [12], or to decelerate and trap molecules 

m 

The reliance on special properties of particular 
molecules has been done away with by the development 
of techniques that combine a static electric field with a 
nonresonant radiative field. The combined fields give 
rise to an amplification effect which occurs for any po- 
lar molecule, as only an anisotropic polarizability, along 
with a permanent dipole moment, is required. This is 
always available in polar molecules. Thus, for a number 
of molecules in their rotational ground state a very weak 
static electric field can convert second-order alignment by 
a laser into a strong first-order orientation that projects 
about 90% of the body-fixed dipole moment on the static 
field direction. If the polar molecule is also paramag- 
netic, combined static electric and magnetic fields yield 
similar amplification effects. So far, the combined-fields 
effects have been worked out for linear polar molecules, 
and corroborated in a number of experiments by the Buck 
[HI IT51 IT5] and Farnik groups [T7] on linear rare-gas clus- 
ters, as well as by the group of Sakai who undertook a 
detailed femtosecond-photoionization study on oriented 

ocs pig ESI . 

Here we analyze the amplification effects that com- 
bined electrostatic and nonresonant radiative fields bring 
out in symmetric top molecules. In our analysis we draw 
on our previous work on linear molecules [50] HI] as well 
as on the work of Kim and Felker [55] , who have treated 
symmetric top molecules in pure nonresonant radiative 
fields. 

For linear molecules, the amplification of the orienta- 
tion arises from the coupling by the electrostatic field of 
the members of quasi-dcgcncratc opposite-parity tunnel- 
ing doublets created by the laser field. Thereby, the per- 
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manent electric dipole interaction creates oriented states 
of indefinite parity [16 . We find that this mechanism 
is also in place for a class of precessing symmetric top 
states in the combined fields. However, another class of 
precessing symmetric top states is prone to undergo a 
sharp orientation via a different mechanism: the laser 
field alone couples states with opposite parity, thereby 
creating indefinite parity states that can then be coupled 
by a weak electrostatic field to the field's direction. 

In Section [TTJ we set the stage by writing down the 
Hamiltonian (in reduced, dimensionless form) for a sym- 
metric top molecule in the combined electrostatic and 
linearly polarized nonresonant radiative fields, and con- 
sider all four symmetry combinations of the polarizability 
and inertia tensors as well as the collinear and perpendic- 
ular (tilted) geometries of the two fields. We derive the 
elements of the Hamiltonian matrix in the symmetric top 
basis set and discuss their symmetry properties. We limit 
ourselves to the regime when the radiative field stays on 
long enough for the system to develop adiabatically. This 
makes it possible to introduce adiabatic labelling of the 
states in the combined fields, sort them out systemati- 
cally, and define their directional characteristics, such as 
orientation and alignment cosines. 

In Section |III| we introduce an effective potential that 
makes it possible to regard the molecular dynamics in 
the combined fields in terms of a 1-D motion. The ef- 
fective potential is an invaluable tool for making sense 
of some of the computational results obtained by solving 
the eigenproblem in question numerically. 

In Section |IV| we present the results proper. First, 
we construct correlation diagrams between the field free 
states of a symmetric top and the harmonic librator, 
which obtains at high fields. Then we turn to the com- 
bined fields, and consider collinear fields and perpendicu- 
lar fields in turn. It is where we discuss the details of the 
two major mechanisms responsible for the amplification 
of the orientation by the combined fields. 

Section [TV] discusses the possibilities of applying the 
combined fields to a swatch of molecules (representing the 
symmetry combinations of the polarizability and inertia 
tensors) and of making use of the orientation achieved in 
a selection of applications. The main conclusions of the 



present work are summarized in Section VI 



II. SYMMETRIC-TOP MOLECULES AND 
THEIR INTERACTIONS WITH STATIC AND 
RADIATIVE ELECTRIC FIELDS 



I are defined such that the principal moments of inertia 
increase in the order I a < I\, < I c . 

The symmetry of the inertia tensor is reflected in the 
symmetry of the polarizability tensor, a, in that the prin- 
cipal axes of the two tensors are collinear. However, a 
prolate or oblate inertia tensor does not necessarily im- 
ply a prolate or oblate polarizability tensor. Four combi- 
nations can be distinguished: (i) I prolate and a prolate, 
i.e., I a < lb — I c and a a < ab — a c with a the figure 
axis; (ii) I prolate and a oblate, i.e., I a < lb = I c and 
a a > otb = a c , with a the figure axis; (iii) I oblate and a 
prolate, i.e., I a = lb < Ic and a a = ab > a c with c the 
figure axis; (iv) I oblate and a oblate, i.e., I a = lb < I c 
and a a = ab < ol c , with c the figure axis. The body-fixed 
permanent electric dipole moment fi of a symmetric-top 
molecule is bound to lie along the molecule's figure axis. 
The symmetry combinations of I and a are summarized 
in Table UJ in terms of the rotational constants 
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and polarizability components parallel, ay 
dicular, aj_, to the figure axis. 



(1) 



and perpen- 



case (i) 


case (ii) 


case (iii) 


case (iv) 


I prolate, 


I prolate, 


I oblate, 


I oblate, 


a prolate 


a prolate 


a. oblate 


a oblate 


I a <I b = I c 


I a <I b = I c 


I a = I b <I c 


I a = I b <Ic 


A>B = C 


A>B = C 


A = B>C 


A=B>C 


a a <a b = a c 


a a >a b = a c 


a a = a b >a c 


a a = a b <a c 


an < aj_ 


aii > a± 


an < aj_ 


ay > aj_ 



TABLE I: Symmetry combinations of the inertia and polariz- 
ability tensors for a polarizable symmetric top molecule. See 
text. 



The rotational Hamiltonian, H r , of a symmetric top 
molecule is given by 



H r = BJ 2 + pBJ 2 z 



(2) 



with J the angular momentum operator, J z its projection 
on the figure axis z (z = a for prolate and z = c for oblate 
top), and 

J (4 - l) > for I prolate , , 

P= \li- 1 )<° for I oblate [6) 



We consider a symmetric-top molecule which is both 
polar and polarizable. The inertia tensor, I, of a sym- 
metric top molecule possesses a three-fold or higher axis 
of rotation symmetry (the figure axis), and is said to 
be prolate or oblate, depending on whether the princi- 
pal moment of inertia about the figure axis is smaller or 
larger than the remaining two principal moments (which 
are equal to one another). The principal axes a, 6, c of 



The symmetric-top molecule is subject to a combina- 
tion of a static electric field, £5, with a nonresonant laser 
field, El- The fields es and El can be tilted with respect 
to one another by an angle, (3. We limit our considera- 
tion to a pulsed plane wave radiation of frequency v and 
time profile g(t) such that 

e|(*) = — 1 9(t) cos 2 (2irvt) (4) 
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where I is the peak laser intensity (in CGS units). We 
assume the oscillation frequency v to be far removed from 
any molecular resonance and much higher than cither t _1 
(with t the pulse duration) or the rotational periods. The 
resulting effective Hamiltonian, H(t), is thus averaged 
over the rapid oscillations. This cancels the interaction 
between [i and 6l (see also ref . [23 , 24 ) and reduces the 
time dependence of El to that of the time profile, 



mum, and eq. (Ill reduces to that for a two-dimensional 



(4(t)) = ^Ig(t) 



Thus the Hamiltonian becomes 



H(t)=H r + V fi + V a (t) 



(5) 



(6) 



where the permanent, V^, and induced, V a , dipole po- 
tentials are given by 



V a (t): 



-Bw cos 9 S (7) 
-BAu(t) cos 2 6 L - Bu>±(t) (8) 



with the dimensionless interaction parameters defined as 
follows 



U) = 

B 

u\\,±_{t) = w\\+g{t) 

2irat\\ ±1 

Wll I = — — 

Be 

Aw = wii — wj_ 
Aoj(t)=aj\i(t) -oj±(t) = Aojg(t) (9e) 



(9a) 
(9b) 

(9c) 
(9d) 



The time-dependent Schrodinger equation correspond- 
ing to Hamiltonian ^ can be cast in a dimensionless 
form 



,4^ = ^(t) 



B dt 



B 



(10) 



which clocks the time in units of h/B, thus defining a 
"short" and a "long" time for any molecule and pulse 
duration. In what follows, we limit our consideration to 
the adiabatic regime, which arises for r ^> h/B. This is 
tantamount to g(t) — > 1, in which case the Hamiltonian 
^ can be written as 



H(t) _ H r 



B 



B 



- uj cos 6s - (Aw cos 2 9 L + u±) (11) 



Hence our task is limited to finding the eigenproperties of 
Hamiltonian (111. For Aw = w = 0, the eigenproperties 



become those of a field-free rotor; the cigenfunctions then 
coincide with the symmetric-top wavef unctions, \ JKM), 
and the eigenvalues become Ejkm/B, with K and M 
the projections of the rotational angular momentum J on 
the figure and space-fixed axis, respectively. In the high- 
field limit, Aw — > ±oo and/or w — > oo, the range of the 
polar angle is confined near the quadratic potential mini- 



angular harmonic oscillator (harmonic librator), see Sub- 
section IIVAI 

If the tilt angle j3 between the field directions is 
nonzero, the relation 

cos 9s — cos (3 cos 9 + sin j3 sin 9 cos tp (12) 

is employed in Hamiltonian ([6]), with 9 = 9l and ip = tpL, 
see Figure [T] 




FIG. 1: Illustration of the angles used in equation |12| for 
substitution for arbitrary field directions. 



If the static and radiative fields are collinear (i.e., 
(3 = or 7r), M and K remain good quantum num- 
bers. Note that except when K — 0, all states are doubly 
degenerate. While the permanent dipole interaction 
mixes states with A J = ±1 (which have opposite par- 
ities), the induced dipole interaction mixes states with 
A J = A ±2 (which have same parities), but, when 
MK ^ 0, also states with AJ — ±1 (which have oppo- 
site parities). Thus either field has the ability to create 
oriented states of mixed parity. 

If the static and radiative fields are not collinear (i.e., 
(3 ^ or 7r), the system no longer possesses cylindrical 
symmetry. The cos <p operator mixes states which differ 
by AM = ±1 and so M ceases to be a good quantum 
number. 

A schematic of the field configurations and molecular 
dipole moments, permanent and induced, is shown is Fig- 
ure H 

The elements of the Hamiltonian matrix in the 
symmetric-top basis set \JKM) possess the following 
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FIG. 2: (color online) Schematic of the configurations of the 
fields and dipoles. An electrostatic, es, and a linearly polar- 
ized radiative, sl, field are considered to be either collinear or 
perpendicular to one another. While the permanent dipole, 
fi, of a symmetric top molecule is always along the figure axis 
(a or c for a prolate or oblate tensor of inertia), the induced 
dipole moment, (at\\ —a±)sL, is directed predominantly along 
the figure axis for an oblate anisotropy of the polarizability 
tensor, an > a±, and perpendicular to it for a prolate polar- 
izability, a 1 1 < a±. See Table [T] and text. 



symmetries 

(J'KM'\H\JKM) = 

(-1) M '- M ( J' —K -M'\H\J -K -M) (13) 
(J' -KM'\H\J -KM) = 

{-1) M '- M {J'K -M'\H\JK -M) (14) 

Since the Hamiltonian has the same diagonal elements 
for the two symmetry representations belonging to +K 
and —K, it follows that they are connected by a unitary 
transformation, U. On the other hand, complex conju- 
gation, JC, of a symmetric top state yields 



K | J KM) = (-1) 



M—K 



J -K -M) 



(15) 



and so we see from eqs. ( 13 1 and ( 14 1 that the +K and 



—K representations are related by the combined opera- 
tion UK., which amounts to time reversal. The two rep- 
resentations have the same eigenenergies and are said to 
be separable-degenerate [25 . We note that the +K and 
—K representations are connected by time-reversal both 
in the absence and presence of an electric field (see also 
ref. [2S]); the time reversal of a given matrix element is 
effected by a multiplication by (— 1) M ~ M . The symme- 
try properties, eqs. (13) and (14 1, are taken advantage 



of when setting up the Hamiltonian matrix. In what fol- 
lows, we concentrate on the case of collinear (i.e., (3 = 
or 7r) and perpendicular fields (i.e., (3 = ir/2). 

We label the states in the field by the good quantum 
number K and the nominal symbols J and M which des- 
ignate the values of the quantum numbers J and M of 
the free-rotor state that adiabatically correlates with the 
state at Aw ^ and/or ui > 0, J -> J and M -> M. We 
condense our notation by taking K to be nonnegative, 



but keep in mind that each state with K ^ is doubly 
degenerate, on account of the +K and —K symmetry 
representations. For collinear fields, M = M . 

In the tilted fields, the adiabatic label of a state de- 
pends on the order in which the parameters u>, Aw and (3 
are turned on, see Subsection |IV C[ As a result, we distin- 



guish among 



J, K,M;uj,Alo,(3 



J,K,M;Au,uj,P 



J, K, M;uj, (3, Alli), and 



J ,K,M; Au>, (3,ujJ states, or 

J, K, M; {p}j for short, with {p} any one of the param- 
eter sets uj, Aoj, (3 or Aid, u, (3 or u, (3, Au or Aw, (3, w. 

The J, K, M;{p}^ states are recognized as coher- 
ent linear superpositions of the field-free symmetric-top 
states, 

J, K,M;{p}) = J2 a JM AI I J > X > M > ( 16 ) 

J,M 

with aj^ M the expansion coefficients. For a given state, 
these depend solely on {p}, 

4^ = aj™ (M) 



(17) 



The orientation and the alignment are given by the 
direction-direction (two-vector) correlation |27j between 
the dipole moment (permanent or induced) and the field 
vector (es or sl)- A direction-direction correlation is 
characterized by a single angle, here by the polar angle Qg 
(for the orientation of the permanent dipole with respect 
to £5) or 6l (for the alignment of the induced dipole with 
respect to el)- The distribution in cither 83 or 8l can 
be described in terms of a series in Legendre polynomials 
and characterized by Legendre moments. The first odd 
Legendre moment of the distribution in 9s is related to 
the expectation value of cos 9s, the orientation cosine, 



(cos9 s )=Y, Y,^Mf*^/ 1 (2J+l) 1 H2J' + l)i 

J'M' JM 

x(-l) M - K 



cos (3 



j i,r 

-KOK 

J 1 J' 

-MOM' 



J 1 J' 
-M -1 M' 



sin/3 



J 1 J 1 

-M 1 M' 



= (j,K,M;{p} 



J,K,M;{p}) (18) 



cos 9 s 



and the first even Legendre moment of the distribution 
in 6l to the expectation value of cos 2 9l, the alignment 
cosine, 
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/„ 2 n \ — „J KM* J KM 

(COS 9 L ) = a J'M a JM 



J'M J 



x(2J' + l)3(-l) 



J,K,M;{p} 



M-K 



J 2J'\f J 2 J' 



-MOM \ -KOK 



COS Ul 



J,K,M;{p}) (19) 



The states with same \K\ and same MK have both the 
same energy, 



E ( J, +K, M; M) = E (j, -K, -M; {p} 



(20) 



and the same directional properties, as follows from the 
Hcllmann-Feynman theorem, 



and 



(cos 9 S ) 



(cos 2 e L ) = - 



m 

dui 



dAoj 



(21) 



(22) 



The angular amplitudes of the permanent and induced 
dipoles are given, respectively, by 



tification algorithm based on a gradual perturbation of 
the field-free symmetric-top states. Instead of relying on 
the eigenvalues, we compare the wavef unctions, which 
arc deemed to belong to the same state when their coef- 
ficients evolve continuously through a crossing as a func- 
tion of the parameters {p}. If \^o) is the state to be 
tracked, one has to calculate its overlap (W^Wo) with 
all the eigenvectors \^k) that pertain to the Hamiltonian 
matrix with the incremented parameters. The state with 
the largest overlap is then taken as the continuation of 
fy . This method has been used for the general problem 
of tilted fields. 



III. EFFECTIVE POTENTIAL 

In order to obtain the most probable spatial distribu- 
tion (geometry), wavefunctions, such as 



J, K, M; M) = *j >K>Mm (0, <p) = 3(0, <p) (25) 



obtained by solving the Schrodinger equation in curvilin- 
ear coordinates, need to be properly spatially weighted 
by a non-unit Jacobian factor, here sin#. Alternatively, 
a wave function, 



|$| 2 = |S| 2 sin( 



#s.o = arccos(cos 9s) 



and 



* [(cos 2 9 S )\ 



(23) 



(24) 



The elements of the Hamiltonian matrix, evaluated in 
the symmetric-top basis set, are listed in Appendix [X] 
The dimension of the Hamiltonian matrix determines the 
accuracy of the eigenproperties computed by diagonaliza- 
tion. For the collinear case, the dimension of the matrix 
is given as J max + 1 — max(|A"|, \M\), with J ranging be- 
tween max (| if |, \M\) and J max - For the states and field 
strengths considered here, a 12 x 12 matrix yields an im- 
provement of less than 0.1% of E/B over a 11 x 11 ma- 
trix. But the convergence depends strongly on the state 
considered. For the perpendicular case, the dimension of 
the matrix is (J ma x + l) 2 — K 2 with J ranging between 
\K\ and J m ax and M between —J and J. Convergence 
within 0.1% can be achieved for the states considered 
with J max = 10. 

The 



J, K, M; {p}j states are not only labeled but also 

identified in our computations by way of their adiabatic 
correlation with the field-free states. The states obtained 
by the diagonalization procedure cannot be identified by 
sorting of the eigenvalues, especially for perpendicular 
fields when the states undergo numerous crossings, gen- 
uine as well as avoided. Therefore, we developed an iden- 



(26) 



with a unit Jacobian can be constructed which gives the 
most probable geometry directly; such a wavefunction is 
an eigenfunction of a Hamiltonian 



H' 

13 



d 2 



U 



(27) 



where U is an effective potential j20J |2T] . Eq. ( 27 1 shows 



that $ corresponds to the solution of a 1-D Schrodinger 
equation for the curvilinear coordinate 9 and for the effec- 
tive potential U . Since $ can only take significant values 
within the range demarcated by U, one can glean the 
geometry from the effective potential and the eigenen- 
ergy. In this way, one can gain insight into the quali- 
tative features of the eigenproblem without the need to 
find the eigenfunctions explicitly. Conversely, one can 
use the concept of the effective potential to organize the 
solutions and to interpret them in geometrical terms. 

For collinear fields, the effective potential takes the 
form 



U-- 



M 2 



sin t 

■U) cos ( 



1 

~ 4 

(Aw cos 



pK 2 ^ — 
2a + ^) 



K 2 - 2M K cos 9 



(28) 



Its first term, symmetric about 9 — n/2, arises due to 
the centrifugal effects and, for \M\ > 0, provides a repul- 
sive contribution competing with the permanent (fourth 
term) and induced (fifth term) interactions. For \K\ > 0, 
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the second term just uniformly shifts the potential, ei- 
ther down when p > (prolate top), or up when p < 
(oblate top). The third term provides a contribution 
which is asymmetric with respect to 9 — tt/2 for precess- 
ing states, i.e., states with MK ^ 0. It is this term which 
is responsible for the first-order Stark effect in symmet- 
ric tops and for the inherent orientation their precessing 
states possess. The fourth term, due to the permanent 
dipole interaction, is asymmetric with respect to 9 = tt/2 
for any state, and accounts for all higher-order Stark ef- 
fects. The fifth, induced-dipolc term, is symmetric about 
9 = tt/2. However, it gives rise to a single well for a pro- 
late (an < a±) and a double- well for a oblate (an > aj_)- 
This is of key importance in determining the energy level 
structure and the directionality of the states bound by 
the wells. 



IV. BEHAVIOR OF THE EIGENSTATES 
A. Correlation diagrams 

In the strong-field limit, a symmetric top molecule 
becomes a harmonic librator whose eigenproperties can 
be obtained in closed form. The eigenenergies in the 
harmonic- librator limit are listed in Tables lU and IIIII and 
used in constructing the correlation diagrams between 
the field- free and the harmonic librator limits, shown in 
Figures [3] (for the permanent dipole interaction, lu — * 00) 
and|4](for the prolate, Aw — > —00, and oblate, Alu — > 00, 
induccd-dipole interaction) . 



= p K 2 + N(2u))i + KM 

+] [3(K - M) 2 -3-N 2 ]-uj 
N=2J-\K + M\ + l 

TABLE II: Eigenenergies for the permanent dipole interaction 
in the harmonic librator limit. See |28l 129) . 

The correlation diagram for the permanent dipole in- 
teraction, Fig. [3J reveals that states with K = split 
into J + 1 doublets, each with the same value of |M|. 
The other states, on the other hand, split into J + \K\ at 
least doubly degenerate states, each of which is charac- 
terized by a value of \K + M\ for a given J. For \K\ < J 
some states are more than doubly degenerate. In the har- 
monic librator limit, the levels are infinitely degenerate 
and are separated by an energy difference of (2a;) 1 ' 2 . 

The correlation diagrams for the induced dipole inter- 
action reveal that states with J < \M\ + \K\ (shown by 
black lines in Fig. [I]) for a oblate and all states for 
a prolate which have same |M_ftT| form degenerate dou- 
blets. In contradistinction, states with J > \M\ + \K\ 
(shown by red and green lines in Fig. [4] for K = and 
K ^ 0, respectively) bound by V a oblate occur as tunnel- 
ing doublets. This behavior reflects a crucial difference 



between the a prolate and a oblate case, namely that 
the induced-dipole potential, V a , is a double well poten- 
tial for a oblate and a single-well potential for a prolate. 

The members of a given tunneling doublet have same 
values of KM and \K\, but J's that differ by ±1. The 
tunneling splitting between the members of a given tun- 
neling doublet decreases with increasing Alu as exp(a — 
frAu; 1 / 2 ), with a,b > 0, rendering a tunneling doublet 
quasi-degenerate at a sufficiently large field strength. In 
the harmonic librator limit, the quasi-degenerate mem- 
bers of a given tunneling doub let coincide with the N + 
and N~ states, see Table III 
with N+ = N- for J < 



Ml 



The N + and N~ states 
+ \K\ always pertain to 



the same J but different KM and so are precluded from 
forming tunneling doublets as they do not interact. In 
the a prolate case, the formation of any tunneling dou- 
blets is barred by the absence of a double well. Note 
that in both the prolate and oblate case, the harmonic 
librator levels are infinitely degenerate and their spacing 
is equal to 2|Aw| 1/2 . 

The correlation diagrams of Figs. [3] and [4] reveal an- 
other key difference between the permanent and induced 
dipole interactions, namely the ordering of levels per- 
taining to the same J. The energies of the levels due 
to Vfj, increase with increasing \K + M\. For V a pro- 
late, they decrease with increasing \M\ for levels with 
\M\ > \K\ while states with \M\ < \K\ have the same 
asymptote. The energy level pattern becomes even more 
complex for V a oblate, which gives rise to \K\ + 1 asymp- 
totes. If \K\ = J, the energy decreases with increasing 
\M\, while for \K\ < J it only decreases for levels with 
\M\ + \K\ > J. All other levels connect alternately to the 
asymptotes with = J or J — 1. This leads to a tan- 
gle of crossings, avoided or not, once the two interactions 
are combined, as will be exemplified below. 



B. Collinear fields 

Figures [5] and [6] display the dependence of the eigenen- 
ergies, panels (a)-(c), orientation cosines, panels (d)-(f), 
and alignment cosines, panels (g)-(i), of the states with 
< J < 3, -1 < MK < 1, and K = 1 on the di- 
mcnsionless parameters lu and Auj that characterize the 
permanent and induced dipole interactions. These states 
were chosen as examples since they well represent the be- 
havior of a symmetric top in the combined fields. The 
two figures show the dependence on Alu for fixed values 
of to and vice versa. Note that negative values of Alu 
correspond to a prolate and positive values to a. oblate. 
The plots were constructed for I prolate with A/B = 2 
but, apart from a constant shift, the curves shown are 
identical with those for I oblate. Thus Figs. [5] and [6] 
represent the entire spectrum of possibilities as classified 
in Table |U 

The left panels of Figs. [5] and [6] show the eigenenergies 
and the orientation and alignment cosines for the cases 
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FIG. 3: Correlation diagram, for the permanent dipole interaction, between the field-free (tj — > 0) symmetric top states | J KM) 
and the harmonic librator states \N) obtained in the high-field limit (uj — ■> oo), see also At intermediate fields, the states 

are labeled by \ JKM 
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FIG. 4: (color online) Correlation diagram, for the induced-dipole interaction, between the field-free (Auj — > 0) symmetric top 
states | J KM) (center) and the harmonic librator states \NKM) obtained in the high-field limit for the prolate, Aui ■— ► — oo (on 
the left), and oblate, Au — > oo (on the right) case. The harmonic librator states are labeled by the librator quantum number 

N and the projection quantum numbers K and M. At intermediate fields, the states are labeled by IjifM^. States that form 

tunneling doublets (only for Auj — > oo) have J > \K\ + \M\ and are shown in color: red for doublets with KM — 0, green for 
doublets with K, M ^ 0. Members of the degenerate doublets have J < \K\ + \M\ and are shown in black. See text. 
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a prolate: Aoj — > — oo 


a oblate: Aoj — * oo 


En >«' m = P K 2 + 

+(2N + 1)(-Aw)s +AP + K 2 - 
-\(2N 2 + 2N + 3)-oj ± 

N = J- |M| for |M| > \K\ 
N = J —\K\ for |M| < \K\ 


+2(N ± + 1)(Acj) L 2 + m 2 2 + K 2 2 ( JV± + 1 ) 2 
— Aoj — ld± — i 

for J < |M| + \K\ 
N ± = 2J — - |M| for KM ^ 

for J > |M| + \K\ 
N± = J - 1 for KM ^ when J —\K + M| is odd 
N ± = J for KM ^ when J - \K + M| is even 

for if or M = 
iV~ = J - 1 when J — M| or J - |K | is odd 
iV+ = J when J - \M\ or J - \K\ is even 



TABLE III: Eigenenergies for the induced dipole interaction in the harmonic librator limit. See also 22 . 



of pure permanent and pure induced dipole interactions, 
respectively. 

For an angle 



7 = arccos 



KM 
J(J + 1) 



(29) 



such that < 7 < 7r/2, the pure permanent dipole inter- 
action, oj > and Auj = 0, panels (a), (d), (g) of Fig. 
[5j produces states whose eigenenergies decrease with in- 
creasing field strength (i.e., the states are high-field seek- 
ing) at all values of oj and their orientation cosines are 
positive, which signifies that the body-fixed dipole mo- 
ment is oriented along £5 (right-way orientation). For 
states with 7r/2 < 7 < tt, the eigenenergies first increase 
with lu (i.e., the states are low-field seeking) and the 
body-fixed dipole is oriented oppositely with respect to 
the direction of the static field ss (the so called wrong- 
way orientation). For states with K = 0, the angle 7 
becomes the tilt angle of the angular momentum vector 
with respect to the field direction, which, for J > 0, is 
given by 



7o 



arccos ■ 



\M\ 



[J(J+1)P 



(30) 



At small oj, states with 70 > 3" 1 / 2 are low-field seeking 
and exhibit the wrong-way orientation, while states with 
70 < arc high-field seeking and right-way oriented. 

At large-enough values of oj, all states, including those 
with K = 0, become high-field seeking and exhibit right- 
way orientation. In any case, the dipole has the lowest 
energy when oriented along the field. Since the asym- 
metric effective potential (28) disfavors angles near 180°, 



and increasingly so with increasing KM, the eigenener- 
gies for a given J decrease with increasing KM. The 
ordering of the levels pertaining to the same J is then 



such that states with higher KM are always lower in 
energy. 

The eigenenergies and orientation and alignment 
cosines for a pure induced- dipole interaction, oj = and 
\Aoj\ > 0, are shown in panels (a), (d), and (g) of Fig. [6] 
The eigenenergies are given by 



E 
B 



B 



= J(J+ 1) + P K 2 - Aoj(cos 2 9) - oj ± (31) 



cf. eqs. (pj) and (111. However, Figs. [6^-c and 12 only 
shows E/B + oj^_ = A, which increase with increasing 
laser intensity for a. prolate and decrease for a oblate. 
Note that both prolate and oblate eigenenergies, E/B = 
A — oj^_, decrease with increasing laser intensity, and so 
all states created by a pure induced-dipole interaction are 
high- field seeking as a result. Note that 



cos' tf = - 



dAu 



(32) 



and thus the alignment cosines are given by the negative 
slopes of the curves shown in Fig. [6] 

In panels (d) and (g) one can see that only the pro- 
cessing states are oriented, and that their orientation 
shows a dependence on the Aoj parameter which qual- 
itatively differs for a oblate and a prolate: for Aoj > 0, 
the orientation is enhanced for states with J < |M| + \ K\ 
and suppressed for states with J > |M| + \K\, while 
for Aw < it tends to be suppressed by the radiative 
field for all states. This behavior follows readily from 
the form of the effective potential, eq. p8| , as shown 
in Figure [T] For the prolate case (Aoj < 0), the po- 
tential becomes increasingly centered at 9 — > 7r/2 with 
increasing field strength, and therefore tends to force the 
body-fixed electric dipole (and thus the figure axis) into 
a direction perpendicular to the field. On the other hand, 
for the oblate case (Aoj > 0), the effective potential pro- 
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FIG. 5: (color online) Dependence, in collinear fields, of the eigenenergies, panels (a)-(c), orientation cosines, panels (d)-(f), and 
alignment cosines, panels (g)-(i), of the states with < J < 3, —1 < MK < 1, and K — 1 on the dimensionless parameter ui 
(which characterizes the permanent dipole interaction with the electrostatic field) for fixed values of the parameter Aui (which 
characterizes the induced-dipole interaction with the radiative field; Acu < for prolate polarizability anisotropy, Alo > for 
oblate polarizability anisotropy). The states are labeled by \ JKM). Note that panels (a), (d), and (g) pertain to the permanent 
dipole interaction alone. See text. 



vides, respectively, a forward (6 — > 0) and a backward 
(6 — ► 7r) well for the N + and N~ states of a tunneling 
doublet (for J > \M\ + \K\) or of a degenerate doublet 
(for J < \M\ + \K\). However, only for the degenerate- 
doublet states does the increasing field strength result in 
an enhanced orientation at £5 = 0. This distinction is 
captured by the effective potential, Fig. [7] whose asym- 
metric forward (for MK > 0) or backward (MK < 0) 
well lures in the J < \M\ + \K\ states. The J > \M\ + \K\ 
states become significantly bound by the V a oblate po- 
tential at Alo values large enough to make them feel the 
double well, which makes the two opposite orientations 
nearly equiprobable. 



Figure [8] shows the dependence of the J-state parity 
mixing on the Aui parameter at ui = (panel a) and 



ui = 10 (panel b). The mixing is captured by a parameter 

e= J2 ("&*) a (33) 

J=2ri+1 

with n an integer. In the absence of even-odd J mixing, 
£ = for J even and £ = 1 for J odd; for a "perfect" 
J-parity mixing, £ — \ for either J even or odd. We see 
that for the sampling of states shown, the non-precessing 
states become parity mixed only when lu > 0. However, 
all precessing states are J-parity mixed as long as Aui ^ 
0. 

The eigenenergies as well as eigenfunctions in the har- 
monic librator limit for both the prolate and oblate po- 
larizability interactions have been found previously by 
Kim and Felker |22j . and we made use of the former in 
constructing the correlation diagram in Fig. [4] We note 
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FIG. 6: (color online) Dependence, in collinear fields, of the eigenenergies, panels (a)-(c), orientation cosines, panels (d)-(f), 
and alignment cosines, panels (g)-(i), of the states with < J < 3, —1 < MK < 1, and K — 1 on the dimensionless parameter 
Aw (which characterizes the induced-dipole interaction with the radiative field; Aw < for prolate polarizability anisotropy, 
Aw > for oblate polarizability anisotropy) for fixed values of the parameter w (which characterizes the permanent dipole 
interaction with the radiative field). The states are labeled by | J KM). Note that panels (a), (d), and (g) pertain to the 
induced-dipole interaction alone. See text. 



that in the prolate case, the eigenenergies and alignment 
cosines, as calculated from Kim and Felker's eigenfunc- 
tions, agree with our numerical calculations for the states 
considered within 4% already at Aw = —50. The pro- 
late harmonic librator eigenfunctions render, however, 
the orientation cosines as equal to zero, which they are 
generally not at any finite value of Aw. 

For the a oblate case, the eigenenergies in the har- 
monic librator limit agree with those obtained numeri- 
cally for the states considered within 5% at Aw = 50. For 
Aw w 350, the difference between the alignment cosines 
obtained numerically and from the analytic eigenfunc- 
tions is less then 3%. We note that for J < \M\ + \K\, 
only one eigenfunction exists, pertaining eithe r to N + for 
KM > or to N~ for KM < 0, cf. Table [TTT] This eigen- 
function is strongly directional, lending a nearly perfect 
right- way orientation to an N + state and a nearly perfect 



wrong- way orientation to an N~ state. These eigenfunc- 
tions pertain to the degenerate doublets. The orientation 
cosines can be obtained from the analytic wavefunctions 
within 0.01% for Aw = 100. For J > \M\ + \K\ (in- 
cluding the cases when either K = or M = 0) both 
the N + and N~ analytic solutions exist, cf. Table |lll| 
and pertain to the tunneling doublets. The degeneracy 
of the doublets in the Aw — > +oo limit precludes using 
the analytic eigenfunctions in calculating the orientation 
cosines. In contrast to the numerical results, the analytic 
solution predicts a strong orientation, which in fact is not 
present, as shown in Fig. [7] The linear combination of 
the analytic eigenfunctions / + and / _ , which is also a 
solution, given as /i, 2 = l/V2(f + ±f-), does not exhibit 
any orientation, just alignment. 

The above behavior of symmetric-top molecules as a 
function of Aw at w = sets the stage for what hap- 
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pens once the static field is turned on and so to > 0. 
For Aoj > 0, the combined electrostatic and radiative 
fields act synergistically, making all states sharply ori- 
ented. For Auj < 0, the orientation either remains nearly 
zero (for states with KM — 0) or tends to vanish (for 
states with KM ^ 0) with increasing |Aw|. 

The synergistic action of the combined fields arises in 
two different ways, depending on whether J < \M\ + \K | 
or J > \M\ + \K\. 

For J > \M\ + \K\, the orientation is due to a cou- 
pling of the members of the tunneling doublets (e.g., the 
|2, 1,1) and |3,1,1) states) by the permanent dipole in- 
teraction. The tunneling doublets occur, and hence this 
mechanism is in place, for Auj > 0. The coupling of 
the tunneling-doublet members by the permanent dipole 
interaction is the more effective the smaller is the level 
splitting between the doublet members (which correlate 
with the N + and N~ in the harmonic librator limit). 
Since the levels of a tunneling doublet can be drawn ar- 
bitrarily close to one another by the induced dipole in- 
teraction, the coupling of its members by even a weak 
electrostatic field can be quite effective, resulting in a 
strong orientation of both states. The wrong- way orien- 
tation of the upper members of the tunneling doublets 
can be converted into a right-way orientation. Such a 
conversion takes place at sufficiently large ui where the 
permanent dipole interaction prevails over the induced 



FIG. 8: (color online) Dependence of the J-parity mixing 
parameter £ on the Auj parameter at uj = (panel a) and 
U! = 10 (panel b). Note that the better the J-parity mixing 
the closer is the £ parameter to 1/2. See text. 



dipole interaction. As noted in previous work |21j . the 
coupling of the tunneling doublets by also arises for 
the non-precessing states, in which case one can speak of 
a pseudo-first-order Stark effect in the combined fields. 
The precessing states show in addition the well known 
first-order Stark effect in the electrostatic field alone, 
which relies on the coupling of states with the same \K\ 
and does not involve any hybridization of J. 

As noted above, the J < |M| + \K\ states are strongly 
oriented by the induced-dipolc interaction alone. Since 
the members of a degenerate doublet that correlate with 
the N + and N~ states (e.g., the |1,1,1) and 1 1, 1, — 1) 
states) have different values of KM, adding an electro- 
static field does not lead to their coupling, as collinear 
Vfj, and V a can only mix states with same KM. However, 
the static field skews the effective potential U, eq. (28 1, 



that enhances the orientation of the right-way oriented 
states (N + ) and, at sufficiently high uj, reverses the ori- 
entation of the initially wrong- way oriented states (N~), 
see below. 

The molecular-axis orientation by which the synergism 
of the static and radiative fields for a oblate manifests 
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itself can be best seen in Fig. [Sf and interpreted with the 
help of the effective potential, Figures [9] and 10 Figs. 
|6|i,e,f provide additional cuts through the same (ui, Aw) 
parameter space. Conspicuously, all states, for a. oblate, 
become right- way oriented at a sufficiently large w, cf. 
Fig. J5f. However, some of the states cither become (e.g., 
|2, 1,0)) or are (e.g., |1, 1, —1)) wrong-way oriented first. 
As u> increases, the |2, 1, —1) state is even seen to become 
right-way oriented, then wrong-way oriented, and finally 
right way oriented again. This behavior is a consequence 
of the different types of coupling that the states in ques- 
tion are subjected to. We'll discuss them in turn. 

The |2, 1,0) state is the upper member of a tunneling 
doublet whose lower member is the |1, 1, 0) state, cf. Fig. 
[4] At u = 0, the |1, 1,0) and |2, 1, 0) states are not ori- 
ented, as is the case for any states with KM = 0. How- 
ever, the value of Aw — 15 is large enough to push the 
levels into quasi-degeneracy, see Fig. [6ji, in which case 
the static field can easily couple them. But at u> <C Aw, 
such a coupling results in localizing the wavefunctions of 
the |1, 1,0) and |2, 1, 0) pair in the forward and backward 
wells, respectively, of the effective potential U, which, for 
w <C Aw, are mainly due to the polarizability interaction. 
As w becomes comparable to Aw, the effective potential 
becomes skewed. The forward well grows deeper at the 
expense of the backward well and the wrong- way oriented 
1 2, 1, 0) state is flushed out into the forward well as a re- 
sult, thus acquiring the right-way orientation. The blue 
effective potentials and wavefunctions in the middle pan- 
els of Figs. [9| and [T0| detail this behavior. We note that 
the lower and upper member of a given tunneling doublet 
is always right- and wrong-way oriented, respectively, at 
lu <C Aw. This is because the coupling by V M makes 
the states to repel each other, whereby the upper level 
is pushed upward and the lower level downward. The 
noted orientation of the two states then immediately fol- 
lows from the Hellmann-Feynman theorem. 

The |1, 1, -1) state has J < \M\ + \K\ and thus, in the 
radiative field alone, is a member of a degenerate dou- 
blet, along with the |1,1,1) state, cf. the blue effective 
potentials and wavefunctions in the upper left and right 
panels of Fig. [9j While the |1, 1, 1) state is always right 
way oriented, the 1 1 , 1 , — 1 ) state is wrong-way oriented 
even at w = thanks to the asymmetry of the effective 
potential due to its angle-dependent third term, propor- 
tional to KM, cf. eq. (28 1. An increase in w removes 



the degeneracy of the doublet and causes the wrong-way 
oriented |1, 1, — 1) state to have a higher energy than the 
right-way oriented |1, 1, 1) state. As w increases, deep- 
ens the forward well and, as a result, the wavefunction 
of the |1, 1, — 1) state rolls over into it, thus making the 
state right-way oriented. 

The |2, 1, —1) state exhibits an even more intricate be- 
havior. Instead of a wrong-way orientation at low w, 
enhanced by the radiative field, the state becomes right 
way oriented first, due to an avoided crossing with the 
1 1 , 1 , — 1} state (whose behavior, sketched above, is, of 
course, also affected by the same avoided crossing). This 



is followed, at increasing w, by a "native" wrong-way ori- 
entation that, at w > 20, is reversed by virtue of the 
deepening forward well of the effective potential, which 
confines the state. 

We note that the tangle of level crossings seen in Figs. 
[5]and[6]that complicate the directional properties of sym- 
metric tops in the combined fields is caused by the re- 
versed ordering of the energy levels due to the permanent 
and induced dipole interactions: as the static and radia- 
tive fields are cranked up, the levels "comb" through one 
another. 



C. Perpendicular fields 

For a tilt angle /3 ^ or tt between the static, £3, and 
radiative, el, fields, the combined-fields problem loses its 
cylindrical symmetry and M ceases to be a good quan- 
tum number. This greatly contributes to the complexity 
of the energy level structure and the directional prop- 
erties of the states produced. At the same time, the 
M-dependent effective potential, eq. (p8|), so useful for 



understanding the directionality of the states produced 
by the collinear fields, cannot be applied to the case of 
perpendicular fields, as M is not defined. 

Figures \TT\ and [T2| show the dependence of the eigenen- 
ergies and of the orientation and alignment cosines on the 
field strength parameters w and Aw for a similar set of 
states as in Figs. [5| and [12| for the collinear fields. 

Figs. [P2fc and d capture well the main patterns of the 
behavior. For a prolate, the interaction with the ra- 
diative field El aligns the body-fixed electric dipole fi 
along the perpendicular static field £s, cf. Fig. [2j For 
a oblate, El aligns fi perpendicular to £5. As a result, 
in perpendicular fields, a prolate yields a strong orien- 
tation whereas a oblate a vanishing one. This is the 
inverse of the situation in collinear fields. However, since 
V a prolate is a single-well potential, the levels lack the 
patterns found for collinear fields for V a oblate. Due to 
the multitude of avoided crossings, the states often switch 
between the right and wrong-way orientation, even over 
tiny ranges of the interaction parameters. Therefore, a 
much finer control of the parameters is needed in the case 
of perpendicular fields in order to preordain a certain 
orientation. For a oblate, the coupling of the different 
states is weak, and the avoided crossings that abound in 
the parallel case, Fig. [6] are almost absent, sec Fig. 12 



The states are essentially all high-field seeking in the 
radiative field for Aw > and low-field seeking for 
Aw < 0. This reflects the repulsive and attractive char- 
acter of the polarizability interaction in the prolate and 
oblate case, respectively. In the oblate case, the states 
shown arc essentially high-field seeking in the static field. 
This means that for Aw > 0, much of the wrong-way ori- 
entation seen, e.g., in Fig. [5}l,f, can be eliminated, see 
Fig. [TT| f. Unfortunately, the elimination of the wrong- 
way orientation happens at the expense of the magnitude 
of the orientation, which remains small. 
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FIG. 9: (color online) Effective potential, U, for K — 1 and MK = — 1 (left panels), Mif = (center panels), and MK = 1 
(right panels) along with the eigenenergies and orientation amplitudes (shown by dots) and squares of the wavefunctions for 
states with J — 1. The column are comprised of panels pertaining to increasing values of u). Red curves correspond to 
Aoj = —15, green curves to Auj = 0, and blue curves to Auj — 15. See text. 



Depending on the relative strength of the two fields, 
the effects of one can dominate those of the other. This 
contrasts with the behavior in the collinear fields where 
even a tiny admixture of the static field can dramatically 
change the behavior of the states due to the radiative 
field (such as the coupling of the tunneling doublets) . 

A detailed comparison of the effect the two fields have 
on a given state is complicated by the dependence of the 
state label on the sequence in which the parameters uj. 
Au> and (3 are varied. This behavior is further analyzed 
below. 

As the dependence of the orientation cosine on the 
u> parameter indicates, see Fig. [TTJ;, the largest ori- 
entation for a prolate polarizability is attained for the 
1 1 , 1,0) state. Other states, such as the |1,1, — 1) state, 
become right-way oriented only for sufficiently large field 
strengths. At u> large, the orientation of all states be- 
comes substantially greater than what is achievable with 
collinear fields, cf. Fig. [Hk For the oblate polarizability, 



the states behave similarly, exhibiting a uniform behav- 



ior. On the other hand, in dependence on Acu, Fig. 12 
the |1, 1, 1) state - instead of the |1, 1, 0) state - exhibits 
the strongest right-way orientation over most of the range 
shown. This is an example of the sequence dependence 
of the state label. The |2, 1,— 1) state shows the most 
dramatic variations with Acu. It has three rather sharp 
turn-around points, where the direction of the dipole mo- 
ment changes. Another example of the sequence depen- 
dence of the label is the |2, 1, 1) state in FigJTTJ;, which 
is wrong-way oriented at u> ~ 10, whereas none of the 
states shown in Figure [T2| i is wrong- way oriented for Auj 
ranging from —20 to 0. 

When Auj 3> u>, quasi-degenerate states are formed, 
similar to the tunneling doublets. However, the electro- 
static field is not able to couple them as a result of which 
the states are only aligned. 

In tilted fields, the label of a given state may depend 
on the sequence in which the fields are switched on and 
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FIG. 10: (color online) Effective potential, U, for if = 1 and MK = -1 (left panels), Mif = (center panels), and MK = 1 
(right panels) along with the eigenenergies and orientation amplitudes (shown by dots) and squares of the wavefunctions for 
states with J — 2. The column are comprised of panels pertaining to increasing values of u). Red curves correspond to 
Auj — —15, green curves to Auj = 0, and blue curves to Auj — 15. See text. 



the tilt angle spanned by the field vectors is varied. Fig- 
ure [13] shows the adiabatic evolution of the states with 
J = 1, K = 1 and M — —1,0,1 for three different se- 
quences leading to a crossing points in Figures [TT] and 
^]atw= 10, Auj = -15 and (3 = tt/2 . The three se- 
quences are: (1) The electrostatic field is turned on to a 
value such that uj — 10; then the laser field is switched on 
to a value such that Aui = —5 • 10~ 5 ; then the tilting of 
the fields is carried out to (3 = tt/2; finally Alj is raised 
in steps of 5 • 10~ 5 up to Auj — — 15. The corresponding 
states are labeled as | J, K, M; uj, (3, Auj). (2) The laser 
field is turned up to a value such that Aw = —15; then 
the electrostatic field is switched on to a value such that 
uj = 5 • 10~ 5 ; the fields are tilted to (3 — tt/2; uj is raised 
(in steps of 5 • 10~ 5 ) to uj = 10. The states are labeled 
as \J,K,M;Auj,(3 7 uj). (3) The laser field is turned up 
to a value such that Auj. Subsequently, the electrostatic 
field is turned on to a value such that uj — 10; finally, the 
fields are tilted to (3 = tt/2. The corresponding states are 



labeled by \ J,K,M;Auj,uj,f3). We note that the states 
| J, K, M; uj, Auj, (3) yield identical results with those ob- 
tained for sequence (3). 
As one can see in Fig. 



13 there are no discontinu- 



ities for any of the calculations that might suggest that 
our state-tracking algorithm has failed. However, the 
three sequences yield different labels for the same state 
at the crossing point. The \\,\,Q;Auj,(3,uj) and the 
[1, 1, 1; uj,(3, Auj) states change their labels for sequence 
2 as compared with sequences 1 and 3. We note that the 
absence of genuine crossings for sequence 1 and 2 pre- 
cludes the possibility that the tracking algorithm jumped 
an avoided crossing. In the middle panel of sequence 2, 
the M — 1 and M — —1 states appear degenerate, but 
they are not. The very weak electric field, necessary to 
give a meaning to the mutual tilting of the fields, intro- 
duces a small splitting. The wavefunctions and thus the 
orientation is different for these two states. The orien- 
tation cosines for all three states of the middle panel for 
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FIG. 11: (color online) Dependence, in perpendicular fields, of the eigenenergies, panels (a)-(b), orientation cosines, panels(c)- 
(d), and alignment cosines, panels (e)-(f), of the states with < J < 3, —1 < MK < 1, and K — 1 on the dimensionless 
parameter ui (that characterizes the permanent dipole interaction with the electrostatic field) for fixed values of the parameter 
Alo (that characterizes the induced-dipole interaction with the radiative field; Auj < for prolate polarizability anisotropy, 
Auj > for oblate polarizability anisotropy). The states are labeled by | J KM). The orientation cosines (cos 9 s) are calculated 
with respect to the electrostatic field and the alignment cosines (cos 2 9l) with respect to the laser field. See text. 



sequence 2 are shown in the panel's inset. 

To our knowledge the above phenomenon of label 
switching has not been previously described. We find 
that it occurs not only for states with the same J and 
different M but also for states with different J. 



Two mechanisms seem likely to be responsible for label 
switching: (i) symmetry breaking; (ii) chaotic behavior 
of the underlying classical system and the concomitant 
singularities of the classical phase space. 

The first mechanism is at hand whenever several inter- 



17 




FIG. 12: (color online) Dependence, in perpendicular fields, of the eigenenergies, panels (a)-(b), orientation cosines, panels(c)- 
(d), and alignment cosines, panels (e)-(f), of the states with < J < 3, —1 < MK < 1, and K — 1 on the dimensionless 
parameter Aui (that characterizes the induced-dipole interaction with the radiative field; Atu < for prolate polarizability 
anisotropy, Alu > for oblate polarizability anisotropy) for fixed values of the parameter uj (that characterizes the permanent- 
dipole interaction with the electrostatic field). The states are labeled by | J KM). The orientation cosines (cos6g) are calculated 
with respect to the electrostatic field and the alignment cosines (cos 2 8l) with respect to the laser field. See text. 



actions are present, only some of which break the sym- 
metry of the system completely. Whereas the symmetry- 
conserving interactions give rise to a block structure of 
the Hamiltonian matrix and so allow for genuine cross- 
ings of states belonging to different blocks, the sym- 



metry breaking interaction leads to just a single block 
and makes all crossings of the states avoided. If the 
symmetry-conserving interactions are turned on first, the 
system may have undergone a genuine crossing before 
the symmetry-breaking interaction is turned on second. 
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FIG. 13: (color online) Comparison of the evolution of states for different path through the parameter space. From top to 
bottom these are: | JKM; u>, f3, Au^, j JKM; Aw, (3, uj^ and j JKM; Auj, u, 0^. In the middle panel the 110 state is interchanged 
with the 11—1 state in the other two panels. We dub this effect 'label switching.' 



However, the opposite is not true, as the symmetry- 
breaking interaction, when turned on first, makes all 
crossings avoided once and for all. 

However, in the case of our system, there seems to be a 
symmetry left, which prevents the system from possess- 
ing avoided crossings only. For sufficiently high densities 
of points in the parameter space, one would expect to 
make any avoided crossings, if present, visible. Never- 
theless, even at the highest point densities, many of the 
crossings have been found to be genuine. For a closed 
path in the parameter space at /3 ^ 0, the evolution of 
states exhibits label switching as well. This is the case 
for closed paths defined in terms of u> & (3 (at fixed Auj) 
as well as in terms of Auj & oj (at fixed /3 - in which case 
the symmetry does not change). Closed paths defined in 
terms of Auj & f3 (at fixed ui) do not lead to label switch- 
ing, which suggests that the permanent dipolc interaction 
plays a major role. From this we infer that label switch- 
ing is not only due to symmetry breaking. We checked 
whether the linear rotor in the combined tilted fields also 
exhibits label switching, and found that it does. 



Arango et al. [31] and Kozin and Roberts [32] have 
shown that a diatomic and a symmetric top in an elec- 
trostatic field alone exhibit classical and quantum mon- 
odromy |33j . In quantum mechanics, monodromy reveals 
itself as a defect in the discrete lattice of states. 

For tilted fields, M ceases to be a good quantum num- 
ber, and its classical analog, m, is no longer a constant 
of motion for a classical diatomic or a symmetric top in 
tilted fields. As a result, the problem becomes nonintc- 
grable. For tilt angles (3 ~ 7r/4, Arango et al. found that 
the diatomic exhibits extensive classical chaos [31], [34] . 
while Kennerly mentions quantum chaos for symmetric 
tops [35j . In order to see monodromy for the tilted fields 
problem, one would wish to construct a quantum lattice. 
However, it is not clear how to go about it in the absence 
of good quantum numbers (apart from K). A possibil- 
ity is to make use of the root-mean-square expectation 
value of M, (M 2 ) 1 / 2 , and construct a quasi-monodromy 
diagram with its aid. Such a diagram can be then com- 
pared with classical results, since (m 2 ) 1 / 2 is well-defined 
in classical calculations [31 . The quasi-monodromy di- 
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agrams that we so constructed strongly depend on the 
ratio of the two field strength parameters. The struc- 
ture of the diagrams for the two types of polarizabil- 
ity anisotropy differs markedly, although it's the elec- 
trostatic field alone that is supposed to determine the 
behavior. The quasi-monodromy diagrams are further 
complicated by the dependence on the representation of 

1/2 

the angular variables, as <^M 2 ) differs in the two rep- 
resentations. Only (M) is well-defined and independent 
of the representation. The quasi-monodromy diagrams 
can be obtain on request. 

Hence, we found that in titled fields switching of the 
adiabatic labels of states for linear molecules as well as 
symmetric tops takes place, where it arises in connec- 
tion with the change of the character of a crossing, from 
genuine to avoided. Monodromy, however, is present for 
both collinear and tilted fields. We note that whether la- 
bel switching could be experimentally observed and uti- 
lized is not clear at this point. 



V. EXAMPLES AND APPLICATIONS 

Table IIVI lists a swatch of molecules that fall under 
the various symmetry combinations of the a and I ten- 
sors, as defined in Table [TJ The table lists the rotational 
constants, dipole moments, polarizability anisotropics as 
well as the values of the interaction parameters uj and 
Auj attained, respectively, at a static field strength of 1 
kV/cm and a laser intensity of 10 12 W/cm 2 . The conver- 
sion factors are also included in the table. While the field 
strength of the electrostatic field of a kV/cm is easy to 
obtain (or sometimes even difficult to avoid), a laser in- 
tensity of a petawatt per cm 2 is somewhat harder to come 
by. However, pulsed laser radiation can be easily focused 
to attain such an intensity, and a nanosecond pulse du- 
ration is generally sufficient to ensure adiabaticity of the 
hybridization process. 

We note that the directional properties displayed in 
Figs. [5] -|T2] should be attainable for most of the molecules 
listed in Table El 

The amplification of molecular orientation in the com- 
bined fields may find a number of new applications. 

In molecule optics [12 , a combination of a pulsed non- 
resonant radiative field with an inhomogcncous electro- 
static field can be expected to give rise to temporally con- 
trolled, state-specific deflections. In an inhomogeneous 
static field produced by an electrode micro-array [38] . 
and a ns laser pulse of 10 12 W/cm 2 , deflections on the 
order of a mrad appear feasible for light molecules. Since 
the synergistic effect of the combined static and radiative 
fields is only in place when both fields are on, the deflec- 
tion of the molecules would be triggered by the presence 
of the ns laser pulse. The sensitivity of the deflection pro- 
cess to the magnitude of the space-fixed electric dipole 
moment and its direction (right- or wrong-way) would 
simultaneously enable state-selection. 

The directional properties of symmetric tops in the 



Molecule 


B 




Aa 


LO 


Auj 




[MHz] 




L J 


@1 

kV/cm 


@10 12 
W/cm 2 


Acetonitrile 


9199 


3.92 


1.89 


0.22 


65.0 


Ammonia 


298500 


1.47 


0.24 


0.0025 


0.24 


Benzene- Ar 


1113 


(0.1) 


-6.1 


(0.05) 


-1735 


Bromomethane 


9568 


1.82 


1.95 


0.10 


64.5 


Chloromethane 


13293 


1.89 


1.69 


0.07 


40.2 


Fluoromcthanc 


10349 


1.85 


0.84 


0.09 


25.7 


Iodomethane 


7501 


1.64 


2.15 


0.11 


90.7 


Trichlormcthane 


3302 


1.04 


-2.68 


0.16 


-257 


Trifluoromcthane 


7501 


1.65 


-0.18 


0.11 


-7.60 



TABLE IV: Values of parameters u> and Auj for choice 
symmetric top molecules whose properties were taken from 
[361 I37| . The conversion factors are: u> = 503.2 fj, [D] 
e s [kV/cnrj/B [MHz] and Auj = 3.1658 x 1CT 7 I [W/cm 2 ] 
Aa [A 3 ]/_B [MHz]. Numbers in parentheses are order-of- 
magnitude estimates. 



combined fields may come handy in the studies of rare- 
gas molecules of the type found by Buck and Farnik. 
Some of the species identified by these authors may in 
fact possess a three-fold or higher axis of rotation sym- 
metry [3"9]. [17] . 

Drawing on the analogy with previous work on spec- 
tral effects in single electrostatic [35], [3D] or radiative 
fields [4], we also note that the combined fields can be 
expected to dramatically modify the spectra of symmet- 
ric top molecules: this is due to a change of both energy 
levels and the transition dipole moments. The latter are 
strongly affected by the directionality of the wavefunc- 
tions, which give rise to their widely varying overlaps. 
However, detailed simulations of such effects in the com- 
bined fields are still wanting. 

Last but not least, molecules in tilted electrostatic and 
radiative fields can serve as prototypical systems for the 
study of quantum chaos. The absence of good quantum 
numbers and the multitude of unstable equilibria suggest 
this possibility. Indeed, as pointed out in the work of the 
Ezra group [31] , [34] . linear molecules in non-collinear 
fields exhibit chaotic dynamics. We hope that our present 
study of symmetric tops in the combined fields will pro- 
vide an impetus for further work on quantum chaos and 
monodromy exhibited by such nonintegrable molecular 
systems. 



VI. CONCLUSIONS 

In our theoretical study of the directional properties of 
symmetric top molecules in combined electrostatic and 
nonresonant radiative fields, we saw that collinear (per- 
pendicular) fields force permanent dipoles of molecules 
with oblate (prolate) polarizability anisotropy into align- 
ment with the static field. We found that the amplifi- 
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cation mechanism that produces highly oriented states 
for linear molecules (K = 0) is also in place for pre- 
cessing states of symmetric tops with J > \K\ + \M\. 
This mechanism is based on the coupling by a collincar 
electrostatic field of the tunneling doublet states created 
by the interaction of the oblate polarizability with a lin- 
early polarized radiative field. The efficacy of this cou- 
pling is enhanced by an increased strength, Auo > 0, of 
the anisotropic polarizability interaction that traps the 
tunneling doublets it creates deeper in a double-well po- 
tential, and draws the members of the doublets closer 
to one another. Apart from this synergistic effect of the 
combined fields, there is another effect in place, but for 
states with J < \K\ + \M\. Such states occur as exactly 
degenerate doublets in the radiative field alone, whose 
both members are strongly but mutually oppositely ori- 
ented along the polarization plane of the field. This ori- 
entation can be manipulated by adding an electrostatic 
field, which easily couples the right-way oriented states 
to its direction, whether this is parallel or perpendicu- 
lar to the radiative field. At a sufficiently large strength, 
to, of the permanent dipole interaction, the wrong-way 
oriented member of the (no longer) degenerate doublet 
becomes also right-way oriented. The above patterns of 
the energy levels are complicated by numerous avoided 
crossings, themselves unavoidable, due to the opposite 
ordering of the energy levels for the permanent and in- 
duced dipole interaction. 

The absence of cylindrical symmetry for perpendicu- 
lar fields is found to preclude the wrong-way orienta- 
tion for the oblate polarizability, causing all states to 
become high-field seeking with respect to the static field. 
The changes of the system's parameters u), Alu, and (5 
cause genuine crossings to become avoided and vice versa. 
This, in turn, causes the eigenstates to follow different 
adiabatic paths through the parameter space and to end 
up with adiabatic labels that depend on the paths taken. 

The amplification of molecular orientation by the 
synergistic action of the combined fields may prove useful 
in molecule optics and in spectroscopy. Monodromy and 
quantum chaos lurk behind the combined fields effects 
whose further analysis may thus be expected to shed 
new light on both. 
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APPENDIX A: APPENDIX 

It is expedient to represent the elements of the Hamil- 
tonian matrix in the basis of the symmetric-top wave- 
functions, | J KM). These are related to the Wigncr func- 



tions D MK via 



\JKM) = {-1) 



M-K 



2J+ 1 
8tt 2 



1/2 



D J _ M _ K (v,0, X )(Al) 



The field operators can be expressed in terms of the 
Legendre polynomials Pj (cos 9) and spherical harmon- 
ics YjM(6,f), which in turn, are related to the Wigncr 
functions by 



D J m {^e,^)=Pj (cosfl) 



1/2 



Y* JM {BM (A2) 
(A3) 



The evaluation of the matrix elements of the Hamil- 
tonian then only requires the application of the "triple 
product over Wigner functions" theorem, 



/ 



Dti 3 KAK)Df} 2K2 (R)Di} iKi (R)dn = 

, J\ J 2 J3 \ I J\ J2 

o7T 



Mi M 2 M 3 I \ K x K 2 K 3 



where R denotes (ip, 0, £). 

(A) For the collinear case, the matrix elements to be 
determined are: 



J'K'M' 



JKM ) = 



(J'K'M' |J(J+ 1) + P K 2 -u±\ J KM) 

-lu(J'K'M'\cos9\ J KM) 

-Alu (J'K'M' |cos 2 6\ JKM) (A4) 



where 



cos 6» = Pi (cos 6>) = D l m (ip,9,£) 



(A5) 



cos 2 9 = i(2P 2 (cos#) + l)= I (2D 2 m (^p, 9, + (A6) 

The properties of the 3j-symbols eliminate interactions 
between states with different M and K. As a result, 
states belonging to certain M and K values can be 
treated in a separate calculation. 



J'KM 



JKM j = 5jj, (J(J + 1) + P K 2 - u ± ) 
-w(2J+l) 1 /2(2J' + l) 1 / 2 



J I J' 

-MOM 

J'=J,J±1 



J I J' 

—K K 



Aw-(2J+l) 1 / 2 (2J' + l) 1 / 



j'=j,j±i 
2 



x(-l) 



M-K 



J 2 J' 
-MOM 



J 2 J' 
-K0K 



.J' = J.J±1,J±2 J' = J,J±,J±2 
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(B) For the case of tilted fields, the matrix elements 
can be expressed either in terms of the static field or the 
laser field. The 'simpler' case is keeping the laser field 
fixed in space and tilting the electric field relative to it. 
Then the cos 9 operator is replaced by: 



cos 9 S — cos (3 cos 9 + sin (3 sin 9 cos ip 



:cos/?£^ + sin/? A /- (£>i 10 



Dlo) 



Different K states do not mix. The matrix elements arc: 



J'KM' 



JKm)=5jj,5 mm > (J(J+1)+ P K 2 -lu ± ) 



■wcos/3(2J+l) 1 / 2 (2J'+l) 1 / 2 

U-K ( J U'\( J U' 

-MOM \-K0K 



j'=j..j±i 



j'=j.j±i 



— lu(2J + 1) 1 / 2 (2J' + l) 1 l 2 (—l) M ~ K sin (3 



J 1 J' 
-M-IM' 



J IJ 1 
—K K 



J' = J,J±1;M'=M+1 J' = JJ±1 



J 1 J' 

-M 1 M' 



J 1J' 
-K K 



J' = J,J±1-M'=M-1 J' = J,7±1 

- Acj(2J+l) 1 ' /2 (2J' + l) 1 / 2 

, u K i J 2J'\( J 2 J' 
-MOM \-K0K 



x (-1)' 



J' = J,J±1,J±2 J' = J,J±,J±2 



If, on the other hand, the laser field is tilted relative 
to the fixed static field, the cos 2 9l operator becomes 

cos 2 8l — cos 2 [3 cos 2 9 + 2 sin (3 cos (3 sin 9 cos 9 cos (p 
+ sin 2 (3 sin 2 9 cos 2 <p 

= 1 + 1 cos 2 pD 2 m + 2^- sm/3cos(3(D 2 _ w -D 2 w ) 



D 2 sin 2 /3 



sin 2 S3 (D 1 
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D, 



20) 



The selection rules then become: J' = J, J ± 1, J ± 2 
and M' = M,M ± 1,M ± 2. The equivalence of the 
two choices of expressing the fields has been numerically 
checked, but only the first, 'simple' choice has been used 
in the calculations. 
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